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experiments. We study the leptonic signatures of such a scenario with tau-sneutrino as 
the NLSP, which is realized in the non-universal Higgs masses scenario. We focus on an 
interesting trilepton signature with two like-sign taus and an electron or a muon of opposite 
sign. The neutralinos and charginos are quite heavy in the model considered, and the 
trilepton signal comes mostly from the slepton-sneutrino production. We identify the relevant 
backgrounds, taking into account tau decays, and devise a set of cuts to optimize this 
trilepton signal. We simulate signal and backgrounds at the LHC with 14 TeV center-of-mass 
energy. Although the sleptons in this model are relatively light, (9(100 GeV), discovery is 
more demanding compared to typical neutralino LSP scenarios. The trilepton signal requires 
large amount of accumulated data, at least ~ 80 fb~^, at the CM energy of 14 TeV. 



1 Introduction 

Supersymmetry (SUSY) is one of the candidates for beyond the standard model (SM) theory 
that is extensively searched for at collider experiments (see e.g. \T]). If supersymmetry does 
exist at the weak scale, some supersymmetric particles are expected to be produced by 
the Large Hadron Collider (LHC) experiment. To discover any supersymmetric signal we 
need a correct theoretical interpretation of the data, and there have been many studies 
on this subject. However, most of these studies assume that the lightest neutralino is the 
stable lightest supersymmetric particle (LSP) u, motivated by its feasibility as a dark matter 
particle [2]. Nonetheless, the neutralino is not the only candidate for dark matter within 
supersymmetry. It has been shown that a gravitino LSP can also be a good candidate for 
dark matter [3]J1]. In this case, due to its very weak interactions the gravitino itself would not 
be seen directly, while the next-to-lightest supersymmetric particle (NLSP) would appear as 
a stable particle at colliders q 

The gravitino dark matter scenario opens up many new phenomenologies with various 
possible NLSPs. If a sneutrino is the NLSP, all other sparticles that are produced at colliders 
would quickly cascade decay to the sneutrino, giving rise to jets and/or leptons along the 
way while the sneutrino itself would yield a large missing energy signature in the detectors. 
Although this is similar to a neutralino LSP scenario, the mass spectrum, production rates 
and branching ratios are in general different, leading to different characteristics. Therefore, 
a dedicated study for the sneutrino NLSP scenario at colliders is justified. Specifically 
we look at a sneutrino NLSP scenario within a supergravity model with non-universal Higgs 
masses. A preliminary study on similar model but within gauge mediated symmetry breaking 
framework has been performed by Covi and Kraml in [5], and recently also analyzed by 
Katz and Tweedie [6] . In this paper, we look into a detailed analysis, involving Monte Carlo 
simulation, of a particular model with tau-sneutrino as the NLSP. We focus on leptonic 
channels in order to find distinguishing signatures of the model, and we study whether the 
signals in such scenarios can be observed at the LHC. 

For hadron colliders, supersymmetric signals can be classified as jets plus missing energy 
(-^t) I2j, jets plus leptons plus missing energy or leptons plus missing energy without a 
jet [8]. Due to the nature of hadron colliders, signals involving jets are expected to have 
higher event rates. However, the standard model backgrounds for these type of events are 
also generally larger. On the other hand, although with relatively small event rates, isolated 



^The LSP is stable if R-parity is conserved, which we also assume in this paper. 

^Due to its long lifetime the NLSP would decay outside of the detectors and appear to be stable. 



multilepton plus missing energy signatures offer relatively clean signals, which in some cases 
can be observed above the SM background. For example, a trilepton signature has been 
proposed as a promising channel to discover supersymmetry with a neutralino as the LSP. 
Motivated by this, and also because the lightest effectively stable particle in our model is 
leptonic, we look at signals with leptons. We found that the trilepton signature in our model 
provides an interesting channel, which can be used to distinguish it from many other models. 
However, search for this signal at the LHC, and hadron colliders in general, is hindered by the 
tau identification problem. On the other hand, inclusive analysis, including jets, shows that 
the LHC at 14 GeV should be capable of discovering the new physics beyond the standard 
model. 

The outline of this paper is as follows. In section 2 we specify our model and its features, 
including the sparticle mass spectrum, decay branching ratios and production rates. In 
section 3 we explore the multilepton signatures of our model. Section 4 consists of discussions 
on the trilepton signals and backgrounds. In section 5 we show the results of our simulation 
analysis for the LHC. We conclude with Section 6. 

2 The Model and Its Features 

For our analysis, we take a specific set of parameters in the Non-Universal Higgs Masses 
(NUHM) model [9]. The free parameters in this model are the universal gaugino mi/2 and 
sfermion masses rriQ, the trilinear coupling Aq (all three at the GUT scale); the ratio of the 
two Higgs vevs tan /3; the Higgs mixing parameter /x; and the CP-odd Higgs mass rriA (at the 
weak scale). It has been shown that sneutrino NLSP is natural in NUHM in the sense that 
it has large parameter space regions allowed by all known constraints from cosmology, dark 
matter and particle physics [10] . Our choice of a model corresponds to NUHM parameters 
tan/3 = 10, mo = 100 GeV, mi/2 = 500 GeV, Aq = 0, /i = 600 GeV and m^ = 2000 GeV. It 
has tau-sneutrino z/^ as the NLSP and relatively light slepton masses, while the squarks and 
gluino are around 1 TeV, still within the reach of the LHC. The full mass spectrum is listed 
in Table [H Note that we use a slightly different value of rrit from the one used in ^U\ . 

We assume that the gravitino mass is lower than rri^^, but we do not need to specify 
its value as it is not relevant here. The lifetime of Uj. depends on the mass gap between 
gravitino G and the tau-sneutrino. However, since the dominant decay channel of tau- 
sneutrino is Ur -^ G + u, the tau-sneutrino would still appear as missing energy even if it 
decays inside the detector. 
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Table 1: The sparticle and Higgs masses of the model we analyze. We assume top pole mass 
rrit = 172.4 GeV [H] and mb{mi,)^^ = 4.25 GeV. The Higgs masses are calculated using 
FeynHiggs [12] . 



Among the lighter sparticles we have the following mass hierarchy 

m^o,m^± > m^o > m^^ > nic^ > m^^ > m^;^ , (1) 

with Ut as the lightest one. Note that the first two generations are mass degenerate in the 
model considered. Throughout this paper we define £ = e, fi, whilst £' = e, /i, r. The decay 
modes for these lighter sparticles are as follows. The chargino can decay as 



,~,± 



xt-^xi + w 



?i + z/£ , iye + 



Ti + Vr, Vr + T 



For the neutralinos the decay modes are 

while for the second lightest neutralino we have the additional decay mode 



(2) 



(3) 



(4) 



although with small branching ratios |j. It is worth noting that the decay modes for xf and 
X2 are similar to that in scenarios with Xi as the LSP. On the other hand, ii exhibits a 
completely different decay pattern, 

h ^ i>£ + f + f , fi+i + T , fi + Ui + U^ , z/^ + £ + z/^ , Dr + l^e + T. (5) 

Note that only 3-body decay channels are open for the selectron/smuon. This is because the 
mass gap between i^ and ui is smaller than rriw and also because of the flavor difference 
between the selectron and the NLSP. The decays in Eq. ([5]) are mediated by virtual W 
{i^ef'f), chargino {uri^^r, fiv^Vr) or neutralino (fi£r, Urii^r) exchange. It is interesting to note 
that the decay mode ii — ^ i't^^^t is highly suppressed because of a destructive interference 
between Xi and X2"^^change contributions. Similarly, for the electron-sneutrino we have 
only 3-body decays 

Ue —^ fi + Ue + T , fi+ i + Ur , Ur + Ue + Ur , Ur + i + T. (6) 

These decays are mediated by virtual chargino (rifz/,-, u^-lr) or neutralino (fii/^r, z/t-z/^z/,-) 
exchange. The decay width of the sneutrino is highly suppressed with respect to the left 
sleptons, see Table [2l Heavier selectron mass provides more phase space and the number of 
accessible decay modes is significantly larger. The stau fi can practically □ decay only to 
the tau sneutrino Ur, 

fi^Ur + f' + f, z>; + z/, + r- , (7) 

where the dominant decay mode is mediated via W [t^rff) and the other one by chargino 
and neutralino. We use SDecay 1.3 [13] to calculate the 2-body decay branching ratios, and 
FeynArts/FormCalc [15] package to calculate the 3-body decay widths U- The branching 
ratios for the decay channels with branching ratios ^ 1% are collected in Table [2l Note that 
the dominant decay mode for z/^ is invisible. 

We calculate the (pair) production rates for the sparticles in our model at the Tevatron 
and at the LHC We assume three center-of-mass (CM) energies for the LHC: 7 TeV, 10 TeV 
and 14 TeV. The results are shown in Table [31 Note that the chargino (xf) and neutralinos 
(Xi 2) are relatively heavy in our model and near the production threshold for the Tevatron. 



•^Note that, X2 can also decay through a loop to Xi + 7 [13,. However, this is subdominant as compared 
to the tree-level two-body decay modes above. 

"^Since the direct decay of stau to gravitino is negligible. 

^The 3-body decays of sleptons and sneutrinos have also been calculated analytically by Kraml and Nhung 
in M- 
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Table 2: Decays and the total widths of Xi , Xi 25 ^Z; ^i- ^'^d fi. Only decays with BR ^1% 
are included. The decay pattern for smuon fi^ is analogous to that for selectron e^. Here 
gu, qd represent m-, c- and d-, s-quarks respectively. Each antiparticle has the same decay 
pattern as its corresponding particle. 



Note, also, that the squarks and gluinos are not produced at the Tevatron because of their 
heavy masses. 



For the light sparticles pair production processes, i.e. Table [3]J^a) and (b), we see that at 
the Tevatron z/^z>* (which is invisible) has the largest cross section due to the light Vr mass, 
followed by r^^v-y and rfz/* (which are the largest visible channels). For the LHC, which is 
a proton-proton collider, t^Vt has the largest cross section, followed by z>^z>* and rf z>*. For 
both colliders gaugino production is subdominant due to their (relatively) heavy masses, and 
in the case of x?, also, by its bino-dominated content. As in most models with neutralino LSP 
(e.g. SPSla [18J) X2Xi^ associated production is the largest among the gauginos, followed by 
XaxX- For comparison, the X2xf production rates for SPSla at the LHC is about 900 fb, 
for 14 TeV CM energy. 

As we can see from the table, squarks and gluinos require large energy because of their 
heavy masses. At 7 TeV, the production of squarks and gluinos is negligibly small. At 
10 TeV, their total production rate is still lower than that of sleptons. At 14 TeV, the gq 
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Table 3: Cross sections in fb for (a) slepton pair, (b) chargino and neutralino pair, and (c) 
squarks and gluino production at the Tevatron and LHC with CM energies 7, 10 and 14 
TeV. The calculation was done with Herwig++ [TTj. Note that squarks and gluino are too 
heavy to be produced at the Tevatron. Here q represents the sum over the light squarks 
u + d+ s + c, while C can be either e or jl. 



becomes important and together with qq provide promising channels for SUSY discovery. 



3 The Leptonic Signatures 



Let us now look at the supersymmetric signals in our model. First, let us focus on the pure 
multilepton plus missing energy signals without associated jetlj. These signals are generated 
from the production of color singlet sparticles, i.e. th^ charginos, neutralinos and sleptons. 
Thus, we look at chargino pair production (x 



xtm B 



neutralino pair production ixlx 



ra^a\ 



lA,J. 



associated chargino-neutralino production {xi X^)i and slepton pair production (^'j^^'^ , C^ i^e', 
i'l^e' and upue', where i' = e,fi,T) as listed in Table [3l From here on, we will implicitly 



^We will consider inclusive searches including jets production in the next section. 
^Same sign chargino pair can only be produced with some associated jets |19] . 



assume the case for the LHC at 14 TeV, unless exphcitly stated otherwise. 

Let us first look closer at the dominant leptonic decay modes for sleptons, sneutrinos, 
the lightest chargino and the second lightest neutralino. The charged sleptons of the first 
two generations can decay directly to Vr as 

il -^ i); + T- + ue, 

it -^ z>, + r+ + z/, . (8) 

Note that the decay channels £^ — ?■ />,- + P,- + £^, i)* + Ur + £^ are suppressed due to the 
cancellations mentioned below Eq. (J5]), in section |2l The selectron/smuon can also decay to 
the respective sneutrino 



-^ Ul + i' + Uii , 



t ^ i>i+ i'^ + i^e' , (9) 



where again i = e, /i and i' = e, fi, r; or with smaller branching ratios to the stau 



-L 



-^ ff +r+ + e^, f + + r- +i^. (10) 



The electron/muon-sneutrino decays mostly invisibly to the tau-sneutrino and neutrinos. 
The largest visible decay mode is 

z/£ — )■ £~ + z/t- + r"*" , 
, ^ £+ + z>; + r-. (11) 



i>* 



They can also decay to stau 



l^E -)► Ti + Ur + 



>; ^ T^ + iy, + i+. (12) 



The leptonic decays of stau are 



f+ ^ 9; + i'+ + z/,, 



T^- ^ Ur+ f ~ + l^i' . (13) 

The second lightest neutralino has much larger production rate (in association with chargino) 
than the lightest neutralino due to its mostly wino content. It decays as 

X2 ^ i'L + ^'^ , ^e + ^e , ^l' + ^i' ■ (14) 



The chargino decays as 



Xt ^ i>e'+i'^, i'^ + iye' 



xr ^ ^; + (' , 4~ + ''i' ■ (15) 

We can classify the pure leptonic signals based on the number of the isolated leptons 
(e, /i and r |j) in the final state as follows: 

A. 1 lepton +-^T- The signals can appear from: 

(1) fft'* (f^Ur) production with the stau decays to z/^ + ue' + ^', where i' could be 
either e, /i or r. 

(2) iJ^Ug {(-t^e) production with the sneutrino decaying invisibly as z/^ — )■ z/,- + z/^ + 1/^, 
while the selectron/smuon decays as £^ — > v'^ + vi + t'. Since the branching ratio 
for selectron/smuon decay to Vr + ^ + Vt is small, the tau final state is dominant. 

(3) xf X2 production with the chargino decays to z>r + r and the neutralino decays to 
Vr + Vt- Again, the tau final state is dominant. 

The standard model backgrounds are coming from direct charged lepton + neutrino 
production through s-channel VT-boson exchange, from single W boson production 
with the cross section of 20 nb [20] , and from WZ with invisible Z with cross section 
3.3 pb |211[22]. These backgrounds are by orders of magnitude larger than the SUSY 
signals which are 0(10 fb). 

B. 2 leptons +-^t (dilepton): The SUSY signals can arise from 

(1) ^^^^ production where each slepton produces one tau through 3-body decay Cl -^ 
z/^ + r + z/£. 

(2) f]^ff production with each stau decaying through 3-body decay mode f i — )■ z/,- + 
z/£/ + £' where d' = e, /x, r. 

(3) z/£Z>| pair production with one of the sneutrino decaying as z>^ — > u^- + i + r 
producing a tau and an electron/muon of opposite signs while the other one 
decays invisibly as z/£ — > z/^- + z/^ + z/,-. 

(4) xrXi" PS'ir production with the charginos decaying to r^z/* and r+z/,- respectively. 



^We will consider tau decays in the next section. 



Contributions from neutralino pair production is suppressed by the small production 
rate. The SM backgrounds come from direct production through 7*,Z* (Drell-Yan); 
from single Z boson production (with cross section of 1.9 nb [2Q]); from ZZ where one 
Z yields a neutrino-antineutrino pair while the other Z yields £'"*"£'" (with cross section 
of 0.3 pb [22l[23]); and from W^W~ production (with cross section of 12.6 pb |2^). 
Again the SM backgrounds are much larger than the SUSY signals. 

C 3 leptons +]^t (trilepton): The SUSY signals can come from 

(1) iJ^Ug {(-t^e.) associated production, followed by dj^ — > u* + u^ + r^ and u} — )■ 
^* _l_ £+ _l_ ^- decays. In this case we have two taus of the same sign and an 
electron or a muon of the opposite sign. 

(2) Xi^X2 associated production, with the chargino decays as xl ~^ ''' + K and 
the neutralino decays as (a) X2 ~^ ^l + ^^ followed by II — > Vr + ^i + t, (b) 
X2 — ^ ^1^ + T^ followed by fi — t- i/^ + z/^/ + £', or (c) X2 ~^ ^i + ^i followed by 

Ug ^ Ur + T + L 

The SM backgrounds for three leptons are from W Z^ and W'y* |j. For the neutralino 
LSP case, in which the dominant channel is through xf X2 associated production, this 
trilepton signature has been studied thouroughly and appears to be a promising channel 
to discover SUSY [26]. In our scenario, however, the xfx2 production is subdominant 
(and certainly insufficient for the Tevatron), and the trilepton signals come mostly 
from d-'ii'i and ^^t'l production. 

It is interesting to notice, however, that for our scenario we can have two taus of same 
sign (i.e. SF+SS), and that the SM background for this is expected to be smaller. In 
this case the SM background receives contribution from three W bosons production 
which has a small cross section. 

D. 4 leptons +IfiT'- The SUSY signals can arise from 

(1) U(,ul production, followed by z>^ — )■ z>^ + r + £ decays. 

(2) (-^(-^ production, with one Cl decaying as £^ — )■ i>^ + d'^ + u^i (where d' = e, /i, r) 
followed by ug — )■ z>T- + £~+r+, while the other slepton decays as £^ — )■ z>T- + r+ + z/£. 

(3) Xi^xr production, with one of the charginos decaying as xf ~^ ^r + '^^ ^^^ ^^^ 
other one decaying through Xi" ^ ^^ + ^'^ followed by z>^ — )■ i/,- + r+ + £^ . 



^At the detector level, there are also some processes that can mimic trilepton signature such as ZZ, tt, 
Drell-Yan and fake leptons |25) . 

9 



The dominant SM background conies from ZZ^ that has cross section of 0.12 pb [23]. 
In our scenario Higgs bosons if and A are quite heavy (~ 2 TeV) and therefore their 
production at the LHC would be suppressed. Moreover, the neutrahnos and charginos 
in our scenario are also relatively heavy. Thus, we do not consider the same kind of 
analysis as done in 



E. 5 leptons +-^t: The SUSY signals can arise from 

(1) £^z/£ {pL^*t) associated production, where the ij^ decays similarly as in the 4- 
leptons case producing 3 leptons while the sneutrino decays as z/£ — )■ z/,- + £ + r. 

(2) Again, neutralino-chargino X2X1 associated production gives subdominant contri- 
bution. Here the neutrahno decays as X2 ^ ^' + ^ followed by £' — )■ i>,- + r + vni 
while the chargino decays as xf — ^ ^'£ + ^ followed by z>^ — )■ z/,- + £ + r. 



'1—1, 



The SM backgrounds are from WZZ [2B], WZY and WYl*- Note that even though 
the SUSY 5-lepton signal has a small rate C(0.1 fb), suppressed by branching ratios 
of ii and ui decays, the SM background is also small. Thus this might also be an 
interesting channel to look at. The question, however, is how much luminosity would 
be needed to receive enough significance. 

4 The Trilepton Signals and Backgrounds 

As mentioned in the previous section, the trilepton signature with a pair of like-sign taus is 
particularly interesting. The signals that we are looking for are T^T^{e, yu)" and T~T~{e, /x)"^, 
which in our SUSY scenario arise mainly from slepton-sneutrino associated production fol- 
lowed by cascade decays (illustrated in Fig. [T]). If the taus and their charges can be identified 
in the detectors then this would provide us with an excellent supersymmetric signal with a 
distinctly larger cross section than the standard model background. In the SM this signature 
can be mimicked, primarily, through the production and decay of three IV-bosons: 

pp -^ W+W+W- -^ T+UrT+uJ-ue , (16) 

and 

pp -^ W-W-W+ -^ T-VrT-Vrtvt , (17) 

respectively. 



"'^"Recall that for this specific model the decay £' -^ Vr + (- + v^ has a very sniaU branching ratio. Having 
this decay channel available, we would have 5-lepton signature with only one tau. 

10 



Vr- 



Vr 


^* 


• 







Figure 1: Example of trilepton signature from slepton-sneutrino associated production. 

In reality, however, taus decay quickly inside the detectors, producing either leptons (i.e. 
e, /i and neutrinos) or jets (plus tau- neutrino). Tau identification could present a problem, 
especially for a hadron collider such as the LHC with high jet multiplicities. If a tau decays 
to e//i then it would be difficult to distinguish it from the electrons/muons produced by 
other processes. On the other hand it is not easy, although not impossible, to identify jets 
that are coming from taus [29||30] . Let us recapitulate on the signals as seen by the detectors: 



(a) eV^(e,/x)^, 

(b) e^e^e^, /i^/x^/i^, 

(c) e^e^/i^, yU^/i^e^, 

(d) r±(e±/i^,/i±e^), 

(e) r±(e±e^,/iV^), 

(f) r±r±(e,/.)^. 

Here r/^ represent a hadronic tau. The ratios are about 

(a) : (6) : (c) : (rf) : (e) : (/) ~ 6 : 3 : 3 : 22 : 22 : 42 (18) 

At this level there is another SM background from the following process: 

PP -^ W+Z -^ T+T+T-U^ (19) 

with the taus decaying either leptonically or hadronically. In addition, there is a background 
from WWW which can produce the leptonic signals directly without going through any tau. 
Signals (a), (b) and (e) also receive backgrounds from WZ and W-^* which produce e/ii 
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directly. Thus, the interesting signals to look at are (c), (d) and (f). Signal (c) provide a 
clear signature, but is suppressed by the branching ratios. Although (d) is quite interesting, 
the signals might be overwhelmed by fake taus. Therefore we concentrate on (f) in our 
analysis, where we look for two tau-jets of same sign and a muon/electron of the opposite 
sign. From here on, we will always mean hadronic tau (r/i) when we say tau (r), unless 
explicitly stated otherwise. 

At the simulation level, we need to consider some additional backgrounds. This is due 
to the fact that there could be some leptons or jets that do not pass the selection criteria, 
resulting in a different signature. For example, we can have ZZ with each Z decaying to a 
pair of taus, then two taus decay hadronically, one tau leptonically while the other one is 
missing. Thus, for t^t^ ji^ signal, we need to include Z-pair, top-pair, and single top - W 
associated production as well. 

The detection of hadronic tau is also not straightforward. Full analysis would require 
detector simulation and tau reconstruction, which are beyond the scope of our paper. To take 
tau identification problem into account we can make an estimate by attaching a detection 
efficiency factor to each hadronic tau, < e/^ < 1. However, to be more precise this factor 
should be taken as a function of transverse momentum and pseudo-rapidity [3ll|32]. Note 
that this factor affects both the signal and the background. For this reason we do not include 
this factor in the histograms shown in the next section below. In addition, the tau charge 
should also be identified correctly for our case. Charge identification is expected to become 
worse for larger tau momentum, although it should not be impossible for the interesting 
range in our model at the LHC [29]. This charge identification can be used to eliminate 
some background events arising from tt, but not entirely. 

At the detector level, there could be additional backgrounds from jets that are misidenti- 
fied as tau's, i.e. fake tau signals. For example, the Wjj which has a much larger production 
rate p3] can be problematic. The rejection rate of fake taus depends on detector's capability 
and is correlated to the tau identification efficiency [29]. If we assume (effective) rejection 
factor of 500, with e^Vejj cross section of 670 pb, we obtain 2.4 fb of fake tau background, 
which is comparable to the SUSY signal. We notice that the missing transverse energy for 
this background is below 200 GeV. 

On the other hand, hadron colliders, such as the LHC, produce many jets in both SUSY 
and SM processes. By looking at trilepton plus any number of jets we would obtain more 
signal events. We start with inclusive search of r^r^fi^ + nj, where n = 0, 1, 2, . . ., and then 
employ cuts to reduce the backgrounds. Our SUSY signal now consists of slepton-sneutrino, 
chargino-neutralino and SUSY QCD. In SUSY QCD, squarks and gluinos are produced and 
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cascade decay to the tau-sneutrino NLSP. We found that this SUSY QCD contribution gives 
large transverse energy to the final states, due to the big gap between the squark sector 
and the slepton sector in our model. Thus this can be used as first evidence of new physics 
beyond the standard model. 

5 Analysis and Results at the LHC 

In this section we study the inclusive trilepton signals at the LHC for a CM energy of 14 TeV. 
The inclusive SUSY signal has been generated with Monte Carlo program Herwig++ 2.4.2 [17] . 
We have included all sparticle pair production processes and all possible 2-body and 3-body 
sparticle decays in the Herwig++ simulations. All background processes have been simulated 
with Monte Carlo program SHERPA 1.2.0 [3l] . For the SHERPA simulations we have used 
COMIX [35] to compute the hard matrix elements. 

After generating events, we apply the following selection criteria: 

1. Jets reconstructed according the anti-Z^-r algorithm with D = 0.7 [36j which are required 
to have 

p^>20GeV, |r/j|<4.5. (20) 

2. A^^ = 1 : Isolated muons with Rf^j > 0.7 . 

3. Nt- = 2 : Isolated like sign taus with R^-j > 0.7 . 

4. The hardest lepton is required to have: 

p^T > 10 GeV , \rie\ < 2.5 . (21) 

5. The two hardest taus in the event are required to have: 

PT^>15GeV, |r/^J<2.5. (22) 

6. Leptons and taus are required to be isolated with 

Re,r, > 0.4 i?,,,,, > 0.4 . (23) 

These form our basic cuts. We have used Rivet 1.2.1 [37] and Fast Jet 2.4.2 [38] in order 
to analyze events according to our prescribed selection criteria. 
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In Fig. [2]we show the transverse momentum distribution for the hardest jet, pip , after the 
basic cuts. It is obvious that the distribution at large px is dominated by contributions from 
SUSY QCD, i.e. from production of squarks and gluinos. Similarly for the second, third and 
fourth hardest jets. This would provide a clear signal of new physics beyond the standard 
model. Note that this feature should also be found for signals with any number of leptons 
in our scenario, and also for other scenarios in which squarks and gluino are much heavier 
than the LSP. Thus, although high px jets indicate new physics, it is not a unique feature 
of our model. 



Pj of the leading jet 



pr of the leading jet 
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Figure 2: The distribution of the hardest jet transverse momentum, p^^ for (a) t^t^ ^ + 
jets and (b) r~T~ ji'^+ jets for SUSY signals and SM background. 



We then apply optimized cuts to enhance the signal-to-background ratio. At this point 
we have two branches of analysis. The main branch is focusing on the leptonic features of 
the SUSY signal, with the following set of cuts: 

A. 1. Veto on 5-jets and more than one jet, i.e. A*";, = and Nj < 1. 

2. Cut on the transverse momentum of the hardest jet in the event above 200 GeV: 
i.e. we require 20 GeV < p^^ < 200 GeV . 

3. Require mmin(/i, r) = min(m(jUi, ri), m(/ii, r2)) < 55 GeV and (pifiji/r) > 1-5 
rads. 

We call this set opt. A for short. It optimizes the SUSY EW signal. For the side branch 
analysis, we have the following set of cuts (opt.B) which is designed to promote the high-p-r 
jets of SUSY QCD signal [39lll0]. 



B. 1. We require A^^ > 2 and ^^ > 200 GeV, and 
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2. A't = E.=lepto„sJets^^ > 300 GeV . 

The results of our simulations are tabulated in Table H] for the i = fi case. The i = e 
case is similar to the fi case, hence is not shown here. Backgrounds for the case of the 
T^T^i~ signature include t^Pt-t^Ut-h^u^ denoted simply as W^W^W~ , t^t^t^Pj. denoted 
as Ziy+, T^T^T^T^ denoted as ZZ, W^[^ T^Vr]W~[-^ fi~Uf,]bbr'^T^ denoted as tiZ, 
t[-^ bT+Ur]i[-^ 6/i~z/^]M/+[^ T+u^] denoted as tiW+, W^[^ r+z/^]iy-[-^ l^'i^f,]bb denoted 
as W~^W~bb, and ^^v^t^Vt-t^t^ denoted as ZW^W~ . Note that the backgrounds from 
top pair production, single top - W boson associated production are already included in 
W^W'bb. We have similar backgrounds for the T~T~i'^, but with the charges conjugated. 

In opt. A, the veto on jets helps to reduce the SM QCD background, in particular ti, 
although it also suppresses the SUSY QCD signal. The cuts on mmin(/i,T) and 0(/i,j/T) are 
used to suppress backgrounds from ZZ and ZW. As we can see from the table, the SUSY 
signal is now comparable to the SM background. In total, it is greater than the background 
for both T^T^fi" and r^r^fi^ cases, but we obtain an improved result for the r'^r^fi^ case. 
With opt.B, on the other hand, SUSY EW signal is suppressed due to the low jet multiplicity. 
We see that after the optimization we obtain a SUSY QCD signal significantly higher than 
the backgrounds. 

We now focus our discussion on the main analysis (i.e. opt. A). In Fig. [3] we show the 
muon transverse momentum distribution pj, after the optimized cuts. The largest back- 
ground comes from ZW. The shape of ZW is following that of SUSY EW, but is is softer. 
However, the signal distribution is larger for smaller pi^, decreasing rapidly with increasing 
Pj.. Therefore there is less incentive to optimize the cut on pj.. 

The transverse momentum distributions of the two taus are shown in Fig. HI Here ti 
is the hardest tau and T2 is the softer tau. The hardest tau momentum peaks at around 
40 GeV and the p^ cut (in the basic cuts) does not reduce the signal much. For the softer 
tau, however, the cut is significant. Increasing the p^ cut from 10 GeV to 30 GeV, for 
example, can reduce the effective cross section by ~ 40%. We choose our cut at 15 GeV, 
which although difficult from experimental point of view should be possible at the LHC. 
Again, we see that it is difficult to reduce the ZW background any further. 

Fig. [5] shows the distribution in the invariant mass of /i and r pair. We note that the 
signal distribution is concentrated at m^^ < 50 GeV. This suggest a bump feature which 
arises from /i and r pair coming from the same decay chain (i.e. from u^, see Eq. ([6])). The 
endpoint of this bump indicates a mass gap of ~ 50 GeV between z/^ and Ur which agrees 
with our mass spectrum. There is, also, a smooth distribution without an endpoint for high 
invariant masses. This arises from pairing the /i with the r that comes from smuon decay. 
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Table 4: Generation characteristics for pp ^ ^ t^t^ + ]pT and pp — )■ /i' 
Tau detection efficiency is not included. 
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Figure 3: The distribution in the (a) /i and (b) /i+ transverse momenta, pj. and pj. 
respectively. Channels giving negligible contribution are not shown here. 
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Figure 4: The distribution in the (a,b) r+ and (c,d) r transverse momenta, p'^ and p^" 
respectively. Channels giving negligible contribution are not shown here. 
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We might be able to cut ZW a little bit more in this case by cutting out around 60 GeV, 
but the gain is not significant. 



invariant mass of 1st T+ and 1st }i 



invariant mass of 2nd T+ and ist /( 



\u 0.04 — 

^ 0-035 h 
^ 0.03 ^ 

J? 0.025 ^ 

5= - 

^ 0.02 ^ 

b : 

■^ 0.015 ^ 

0.01 H 

0.005 ^ 

o E- 



n 



— l+v 
-- ZW+ 

• ■• zz 




150 200 

'"tI,;,! [GeV] 




200 

,,1 [GeV] 



(b) 



invariant mass of ist t and ist ^ 




150 200 

niri,;,! [GeV] 



150 200 

mT2,;,i [GeV] 



(d) 



Figure 5: The distribution in the (a,b) fi'r'^ and (c,d) /i^r^ invariant masses. Here ri is 
the hardest tau, and T2 is the second hardest. Channels giving negligible contribution are 
not shown here. 



In Fig. |6]we show the invariant mass distribution of tt pair. There is no endpoint feature 
seen in this plot, suggesting that the two taus always come from the opposite decay chains. 
On the other hand, notice that the invariant mass distribution peaks at around 50 GeV, 
indicating that both taus are coming from decays of weak scale particles. 

In all of the plots above, the tau efficiency factor e has not been included. Note that 
since both signal and background are affected by e, including this factor would only rescale 
the distribution height but would not change the ratio between signal and background lj. 



"'^"'^Note, Jiowever, that e varies with respect to some observables such as the tau transverse momentum pip 
and rapidity. Therefore the rescaUng is not constant. 
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Figure 6: The distribution in the (a) t^t~^ and (b) r r invariant masses. Channels giving 
neghgible contribution are not shown here. 



except for the fake tau rate which is not included in the plots. Here we assume that the fake 
tau rate, which depends on the real data analysis, can be kept under control. 
Significance can be estimated as follows 

S 



<ys 



^/sTb 



■ Ceff ■ 




(24) 



1/(75 + aB 

where eetr is the effective tau identification efficiency factor over the whole spectrum and / C 
is an integrated luminosity. Here we have used the fact that there are two taus in our signal, 
and assumed (for simplicity) the same effective efficiency factor e^s for both taus. The tau 
charge identification efficiency is implicitly included in egg, i.e. e = e^Ccharge- 

Recalling the effective cross sections for the T^Th^, signal after the cuts as summarized 
in Table HI we find that the total effective cross sections are cropt.A(SUSY) ~ 3.2 fb and 
'7opt.A(SM) ~ 1.9 fb respectively, including both t^t^ jj,~ and t^t^ ji'^ . We see that for 5-0" 
discovery level, the required integrated luminosity is 

>C(5a)opt.A ^ 12.5/e,\ (fb"^) . (25) 

Taking eeg = 0.4 [29], for example, we find that 5-0" discovery requires about 80 fb~^ of data. 
On the other hand, for p^ analysis (opt.B), we have o"opt.B(SUSY) ~ 7.6 fb and o"opt.B(SM) c 
0.7 fb respectively, leading to 

y'£(5a)opt.B^3.6/e,\(fb-i). (26) 

Thus, with eeff = 0.4, 5-0" level in this case requires only about 23 fb~^ of data. Note that 
this is not necessarily the most promising channel to look for high-p^ jets, i.e. we should 
also compare it with pure jets channel etc, each of which requires a separate set of analysis. 
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Here, we have assumed 14 TeV CM energy for our analysis. Even using the highest 
energy expected at the LHC we see that we need a significant amount of data for discovery. 
We can deduce from Table [3] that it would be very difficult with 10 TeV CM energy and 
practically impossible with 7 TeV. Thus, we hope that the LHC can overcome its technical 
difficulties and reach the original designed energy of 14 TeV. 

6 Conclusion 

We have studied the leptonic signatures of a model in which tau-sneutrino is the effectively 
stable lightest supersymmetric particle at the LHC. The model that we consider has relatively 
heavy charginos and neutralinos, and relatively light sleptons and sneutrinos. The cross 
sections for pure leptonic signals are generally small, partly due to the fact that neutralino- 
chargino associated production in this model is suppressed by the heavy gaugino masses. 
Nevertheless, we find that the trilepton signature is still interesting to look at. It consists of 
a signal with two like sign taus and one electron or muon of opposite sign, coming from ClVi 
production. 

We employ an inclusive search strategy in generating signals and use a set of cuts to look 
at this particular signature. At the Tevatron, the sparticle production rates are too small to 
yield any observable supersymmetric signal in our scenario. At the LHC, sufficiently large 
CM energy is still required. With 14 TeV, and the optimized cuts, we can obtain 5-cr SUSY 
trilepton signals after ~ 80 fb^^ integrated luminosity. We, also, investigated the leptons + 
jets signatures, and noticed that we can use pj- cut to observe new physics signal above the 
standard model background. In our case, for pj, ^ 200 GeV, where ji is the hardest jet, 
SUSY QCD is dominant due to the large mass gap between the squark sector and the slepton 
sector. In this way we can obtain a significant signal-to-background ratio after 23 fb~^ of 
integrated luminosity. However, this does not tell us much about the underlying model, since 
this is generally true for any supersymmetric model with heavy squarks and gluino. Note 
that these are optimistic estimates, assuming that we can suppress the fake-tau event rate. 

Our study suggests that the search for supersymmetry can be quite challenging, depend- 
ing on the specific supersymmetric model. This is especially true when we want to look 
into the detailed characteristics of the model. Even though the slepton spectrum is rela- 
tively light, around 100 GeV, we still need large amount of data and high energy to see a 
significant excess of signal over background. 

If this scenario is realized in nature, a big challenge in the data analysis at the LHC 
would come from tau reconstruction and identification, as well as rejection of fake taus. 
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As tau can appear copiously in many models beyond SM, we might need new methods in 
this aspect. Indeed, there are ongoing efforts to alleviate these problems [H]. Nevertheless, 
hadron colliders, such as the LHC, might not be sufficient to explore the physics beyond 
the Standard Model. In this case, a lepton collider such as the proposed e~^e~ International 
Linear Collider (ILC) [12] would help. For the ILC, the tau signal would be much cleaner. 
Only then it would be possible to probe the model further, for example by reconstructing 
masses and measuring couplings. 
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